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1. Introduction. Center manifolds are an important tool for the local analysis 
of dynamical systems. In this paper we develop a methodology to prove the existence 
of center manifolds in a "large" neighbourhood of the equilibrium point. The method 
involves the use of normal forms, topological results, and computer assisted compu- 
tations. The novelty of our approach is that it provides explicit rigorous bounds on 
the size and location of the manifold for a given dynamical system. Moreover, under 
appropriate hypothesis we prove that the manifold is unique. 

In contrast, the classical center manifold theorems show existence of a manifold 
in some neighbourhood, but they do not readily provide information on the size of 
this neighbourhood. Also, the classical normal form theorems construct an accurate 
approximation to the dynamics in a neighbourhood, but the normal form is usually 
not convergent. Sometimes the normal form does converge, but we lack information 
on its domain of convergence. 

To show the power of our methodology, in the second part of this paper we 
prove existence and uniqueness of the center manifold in a practical domain around 
an equilibrium point of the celebrated Restricted Three Body Problem (RTBP). By 
practical we mean that such domain possibly could be used for realistic space mission 
design, since it is not too small. To our knowledge, this is the first proof of existence 
of the center manifold in a practical domain for the RTBP. 

For the rest of this introduction, we define the center manifold and mention some 
previous results related to this paper. Finally we explain how this paper is organized 
into sections. 

Definition 1.1. Consider a differential equation on W 1 

x = Ax + f(x), (1.1) 

where A is linear and f has no constant or linear terms. The origin is a fixed point. 
Let R n = E c ® E u E s be the usual decomposition into the center, unstable, and 
stable invariant subspaces with respect to A. 
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A center manifold W c is an invariant manifold of the flow of (1.1), tangent to 
E c at the origin, and of the form 

W c = {(9, X (e)): 9eU}, 

where x '■ E c — >■ E u E s is a C k function, and U is an open neighborhood of in 
E c . 

We are naturally lead to study the flow in the center manifold. The center mani- 
fold approach has the advantage that this reduced problem is a dynamical system on 
a lower-dimensional manifold (of the same dimension as E c ). The reduced problem 
contains crucial information of the full problem (1.1). The qualitative behavior of the 
flow on the center manifold completely determines the behavior of the full flow around 
the fixed point [Car]. Also, every center manifold contains all globally bounded so- 
lutions (e.g. fixed points and periodic orbits) which are close enough to the origin 
[Sij]. 

Let us now mention some results related to this paper. The existence of center 
manifolds is discussed in many dynamical systems books, for instance in [GH, CH], 
and the monograph [Car]. The subtle properties of center manifolds such as (non)- 
uniqueness, tangency, limited differentiability, and (non)-analyticity are discussed 
in [Sij]. 

Lyapunov [ ] studied the case in which the linear operator A in equation (1.1) has 
a simple pair of eigenvalues ±wi. He proved the existence of a center manifold filled 
with an analytic one-parameter family of periodic orbits. The main hypothesis are 
the presence of an integral of motion and a nonresonance condition. Such situation 
arises for the equilibrium point of the Restricted Three Body Problem that we study 
in the second part of this paper. Lyapunov's theorem applies to the Restricted Three 
Body Problem (cf. [SM] §18), but again is only local and does not readily provide 
estimates on its domain of validity. 

Normal forms make a very powerful and general technique to approximate local 
dynamics, including the center manifold, and stable/unstable manifolds of a fixed 
point. It is also a classical subject discussed in many dynamical system books, for 
instance [Mil] (Hamiltonian systems), [CH] (general differential equations), [GH], and 
the monograph [Mil]. 

Given their usefulness, normal forms have been applied to approximate the center 
manifold around the equilibrium points of the planar Restricted Three Body Problem 
[CM, CDMR], and the spatial RTBP [JM, DMR]. In particular, our implementation 
of normal forms is based on [J]. This technique has important applications in space 
mission design [GJSM, GKLMMR] and diffusion estimates [JS, JV]. 

Regarding the planar RTBP, we would also like to mention the numerical explo- 
rations of Broucke [B] , where he performed an extensive study of different families of 
periodic orbits. In particular, he finds a family of numerical periodic orbits around 
the same equilibrium point that we study in this paper. The family extends up to 
a very large neighborhood of the equilibrium point (much larger than our rigorous 
result). 

The paper is organized as follows. In Section 2 we give the setup of the problem 
and state our main theorem (Theorem 2.4). Assumptions of the theorem are based on 
estimates on the derivatives of the vector field within the investigated region. Based 
on these the existence of an invariant manifold is established. In Section 3 we give a 
topological proof of the existence of an invariant manifold for maps with saddle-center- 
type properties. In Section 4 we use the result obtained for maps to prove Theorem 



CENTER MANIFOLDS IN A PRACTICAL DOMAIN 



3 



2.4. In Section 5 we apply our Theorem 2.4 to prove the existence of a center manifold 
around an equilibrium point Li in the RTBP. To do so we first introduce the problem 
and present a procedure of transforming the system into a normal form. We then 
discuss how normal forms provide very accurate approximations of center manifolds. 
Finally we combine Theorem 2.4 and normal forms with rigorous interval arithmetic 
based computer assisted computations to prove the existence of the manifold. Section 
6 contains concluding remarks and an outline of future work. 

2. Setup. We will consider the following problem. Let F : E n — > E™ and 

x' = F(x) (2.1) 

be an ODE (we impose the usual assumptions implying existence and uniqueness of 
solutions) with a fixed point Xo and an integral of motion H : R™ — > E. By this we 
mean that for any solution q(t) of (2.1) we have 

H(q(t)) = c, (2.2) 

where c is some constants dependent on the initial condition q(0). Since in our appli- 
cations we shall deal with the restricted three body problem, which is a Hamiltonian 
system where H is the Hamiltonian, we shall refer to H as the energy from now on. 
We shall use a notation $(i,x) for the flow induced by (2.1). 

2.1. Well aligned coordinates . We will investigate the dynamics of (2.1) in 
some compact set D, contained in an open subset U of E n , such that the fixed point 
xo G D, and whose image by a diffeomorphism 

<t> : U -> 4>{U) C E™ (2.3) 

is 

4>(D) =Di = B? xB r u x Bl, (2.4) 

where B\ (for i £ {c, u, s}) stand for z-dimensional closed balls around zero of radius 
r. We assume that n = c + u + s. We will refer to p = 0(x) as the aligned coordinates. 
In these coordinates we will use a notation p = (9,x,y) with 9 g B^, x £ B T U and 
y G B r s . We will refer to 9 as the central coordinate, to x as the unstable coordinate 
and to y as the stable coordinate (the subscripts c, u, s standing for central, unstable 
and stable respectively). 

The motivation behind the above setup is the following. We will search for a 
center manifold of (2.1) homeomorphic to a c-dimensional disc inside the set D. Such 
manifolds have associated stable and unstable vector bundles, which in the coordi- 
nate system <j) are given approximately by the coordinates of the balls B r s and B r u 
respectively. We do not assume though that the coordinates x and y align exactly 
with directions of hyperbolic expansion and contraction. It will turn out that it is 
enough that they point roughly in these directions. The remaining coordinates 9 are 
the central coordinates of our system. We need to have a good guess on where the 
center manifold is. This guess is given by x {0}) C E™. The change of co- 

ordinates <j) can be obtained from some non-rigorous numerical computation (in our 
application for the RTBP - normal forms). It is important to emphasise that we will 
not assume that (f>~ 1 (B^ x {0}) is invariant under the flow (2.1). Allowing for errors, 
we expect the true manifold to lie in _1 (i3 ( f- x B r u x B r s ). This means that we take 
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an enclosure of radius r of our initial guess and look for the invariant manifold in this 
enclosure. 

We will search for the part of the center manifold with energy H < h for some 
h G M. We assume that the center coordinate is well aligned with the energy H in 
the sense that we have 



Htf-HB?- V x K x BD) <h<H {cf>-\dB? x B r u x B r s )) , 



(2.5) 



for some v > (here we use a notation dA to denote the boundary of a set A). 

Our detection of the center manifold in the RTBP is going to be carried out in 
two stages. First we obtain as a change of coordinates into a normal form, after 
which wc shall employ our topological theorem (Theorem 2.4) to prove the existence 
of the manifold. 

2.2. Local bounds on the vector field and the statement of the main 
result. We are now ready to state the main assumptions needed for our method. 
These will be expressed in terms of local bounds on the derivative of the vector 
field (2.1). First let us introduce a notation for the vector field in the aligned 
coordinates i.e. 

F^(p) = D^- 1 (p))F^- 1 (p)), (2.6) 
and a notation [dF^(N)] for an interval enclosure of the derivative on a set N C 



[dF^{N)\ 



Ae 



\An G 



inf 



dFf 



dFf 



.... , -,sup— — 

p£N dpj peA r dpj 



for all i,j = 1, 



For any point p = (8, 0, 0) from B^ x {0} x {0} we define a set 
N p :=B c (9,p) xB r u x B r s n D h 



(2.7) 



where B c (8,p) is a c dimensional ball of radius p > centred at 9. We introduce the 
following notations for the bound on the derivatives of F^ on the sets N p 



[dF*(N p )] C 




(2.8) 



Here A, B, C,e c , e m , e s and e u are interval matrices, that is matrices with interval 
coefficients. Here we slightly abuse notations since the pairs of matrices e c and e m 
need not be equal; they even have different dimension when u ^ s. We use the same 
notation since later on we shall assume uniform bounds for both of matrices e c and 
both e m . Let us also note that the bounds A, B, C,e c , e m , e s and e u may be different 
for different p. We do not indicate this in our notations to keep them relatively simple. 

Remark 2.1. If the system possesses a center manifold and the adjusted coor- 
dinates are well aligned in the sense of section 2.1, then the interval matrices £j in 
(2.8), with i G {c, m, s,u} should turn out to be small. The matrices A,B,C are 
the bounds on derivatives of the vector field in the unstable, stable and central direc- 
tions respectively. If the alignment of our coordinates is correct then we expect the 
contraction/ expansion rates associated with C to be weaker than for A and B. 

We will use the following notations to express our assumptions about [dF^(N p )]. 
Let S u , 5 s , c u , c s , Ei > denote contraction/expansion rates, such that for any matrix 
A G A, B G B, ei G e, for i G {m, c, u, s}, we have 
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M{x T Ax : \\x\\ = 1} > S u , (2.9) 

S up{y T By:\\y\\ = l}<-S s , (2.10) 

c s < mi{0 T C9 : \\6\\ = 1} < sup{6> T C6> : \\9\\ = 1} < c", (2.11) 

||e,||<£j for i G {m, c, u, s}. (2-12) 

Once again, c s , c u , /x, 5" and <5 S can depend on p. 
Let j,ah,otv, Ph, ftv > be constants such that 

a/, > a„ and /?„ > /?/,, (2-13) 

and such that the radius p considered for the central part of the sets N p satisfies 

P>r x f^, p>rM, (2.14) 
V 7 V 7 

where r is the radius of the balls B T U and B r s in (2.4). Let us define the following- 
constants 

2 V 7 7 / 

„forw ._ X u 1 ( „ , _ , 7 _ , 0h 



^ - o £ ™ + e " + — e = + — £ * . ( 2 - 15 ) 

^forw . = _ja + 1 ( + 0* + 1 + £ 
* \ Ph Ph 



s I I 



K c . C ( £ m ~t- £ m -\- 2s ( 

2 V 7 7 

K back , = X u _ 1 ^ + £u + + A, £s | (2 16) 

2 V a f a« ' 



2 V & A; 

The superscripts "forw" and "back" in the above constants come from the fact that 
they shall be associated with estimates on the dynamics induced by the vector field 
(2.6), for forward and backward evolution in time respectively. At this stage the 
subscripts v and h in constants a and (3 do not have an intuitive meaning. During 
the course of the proof they shall be associated with horizontal and vertical slopes of 
constructed invariant manifolds (hence h for "horizontal" and v for "vertical"), and 
then their meaning will become more natural. 

Remark 2.2. Even though coefficients (2.15), (2.16) are technical in nature, 
they have a quite natural interpretation in terms of the dynamics of the system. 
The estimates k % c ,k % u ,k\ for i € {forw, back} are essentially estimates on the con- 
traction/expansion rates associated with the center, unstable and stable coordinates 
respectively. These estimates take into account errors Si for i G {s,u,c, m} in the 
setup of coordinates. Note that when our coordinates are perfectly aligned with the 
dynamics, then Si = for i € {s, u, c, m}, and in turn 



.forw back rs back s /ora u ^forw back r)i 

— K s — O , K c — C , tv c — C , K u — h u — , 
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which are the bounds on the derivative of the vector field in the unstable, stable and 
center directions given in (2.9), (2.10), (2.11). The key assumptions of Theorem 2.4 
are (2.17) and (2.18). In particular, these assumptions imply 

back _^ back „forw ~ /odd 

which is equivalent to assuming that the dynamics in the center coordinate is weaker 
than dynamics in the stable and unstable directions. These are classical assumptions 
for center manifold theorems (See [GHJ, for instance). 

Remark 2.3. We have certain freedom of choice for the constants 7, a^, a v , 
j3h, (3 V - They offer flexibility when verifying assumptions of Theorem 2.4- During the 
course of the proof of Theorem 2.4 it will turn out that they also give Lipschitz type 
bounds L c = mMah JZ,p^ j' L * = ^ m ax( 7 ,^), L u = ^max( 7 ,a„) for 
our center, stable and unstable manifolds respectively (for more details see Corollary 

4.4). 

We are now ready to state our main tool for detection of center manifolds. 

Theorem 2.4. (Main Theorem) Let del. Assume that (2.5) holds for some 
v > 0. Assume also that for any p e B^ x {0} x {0}, for the constants n{° rw , 
4 orw , Kc ack , K b u ack , Ks ack , e u , e s , S u , 6 s computed on a set N p (defined by (2.7)) the 
following inequalities hold: 

K f s orw < K { orw , < (2.17) 

K b s aCk < 0, K b s aCk < K b c ack , K b u ack , (2.18) 

and also that there exist E u , E s > such that for any q £ N p n (B^ x {0} x {0}) 

||7r B i^(?)|| < rE u , \\TT y F*{q)\\ < rE s , (2.19) 

and 

E u +e u <S u , (2.20) 
E s + e s <S s . (2.21) 

If above assumptions hold, then there exists a C° function 

X : B?- v -> D$ 

such that 

1. For any 9 G B^~ v we have TTgx(&) = 9 o,nd 

S D forallteR. 

2. If for some x e c/)^ 1 (B^ v xBJx B r s ) we have 

<£>(t,x)eD forallteR 

then there exists a 9 € B^~ v such that x = (j>^ 1 (x(^))- 

In subsequent sections we shall present a proof of this theorem building up aux- 
iliary results along the way. Before we move on to these results let us make a couple 
of comments on the result. 

Remark 2.5. Theorem 2.4 establishes uniqueness of the invariant manifold. This 
is not a typical scenario in case of center manifolds which are usually not unique. 
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Uniqueness in our case follows from condition (2.5), which by our construction will 
ensure that for any point from our center manifold a trajectory starting from it cannot 
leave the set D. This means that dynamics on the center manifold with H < h is 
contained in a compact set. This is the underlying reason that allows us to obtain 
uniqueness. 

Remark 2.6. The main strength of our result lies in the fact that it allows us to 
easily obtain explicit bounds for the position and size of the manifold. The center man- 
ifold is contained in D — <f>~ l {B^ x B r u x B r s ). Since the manifold is a graph ofx> from 
point 1. of Theorem 2.4 we know that it is of the form (f)" 1 {(9,Tr x y x(0)) \8 £ B^~ v } , 
which ensures that it "fills in" the set D nontrivially. In contrast, the classical center 
manifold theorem does not provide such explicit bounds. 

Remark 2.7. In principle, one could derive some explicit analytic bounds using 
e.g. the "method of majorands" explained in the book of Siegel-Moser [SM]. However, 
to guarantee existence of the center manifold in a neighborhood of the equilibrium point 
that is not too small, one would require a substantial amount of very careful estimates. 

Remark 2.8. It is important to remark that our result only establishes continuity 
(together with Lipschitz type conditions) of the center manifold. The center manifold 
theorem clearly indicates that in a sufficiently small neighbourhoods of a saddle-center 
fixed point we should have higher order smoothness. We believe though that similar 
in spirit assumptions to those of Theorem 2.4 should imply higher order smoothness. 
This will be the subject of forthcoming work. The result obtained so far should be 
regarded as a first step towards this end. 

In our application for the RTBP, in a neighbourhood sufficiently close to the equi- 
librium point, our manifold shall inherit all regularity which follows from the center 
manifold theorem (see Remark 5.1). 

Let us finish the section with a final comment. In order to verify assumptions of 
Theorem 2.4 it is sufficient to consider some finite covering {IJie/ Ui} of the set 
and to verify bounds on local derivatives on sets Ui. It is not necessary to consider 
an infinite number of points p and their associated sets N p , as long as for any p £ 
Bf x {0} x {0} we have N p C U { for some i £ I. This makes assumptions of Theorem 
2.4 verifiable in practice using rigorous computer assisted tools. 

3. Topological approach to centre manifolds for maps. In this section we 
will state some preliminary results, which will next be used for the proof of Theorem 
2.4 in Section 4. The results will be stated for maps instead of flows. In Section 4 
we will take a time shift along a trajectory map for the flow generated by (2.1) and 
apply the results to it. The main result of this section is Theorem 3.7. The result is 
in the spirit of versions of normally hyperbolic invariant manifold theorems obtained 
in [Ca], [( 5] and [CS]. The main difference is that we do not deal with a normally 
hyperbolic manifold without boundary, but with a selected part of a centre manifold 
(homeomorphic to a disc) with a boundary. In this section the fact that the dynamics 
does not diffuse through the boundary along the centre coordinate is imposed by 
assumption. This assumption will later follow from assuming that (2.2), (2.5) hold 
for (2.1). 

We now give the setup for maps. Let D C U C K n , the change of coordinates 
4> : U —> 4>(U), and = <j>(D), be as in Section 2.1. We consider a dynamical system 
given by a smooth invertible map /:[/—> U. In adjusted coordinates we denote the 
map as := <f> o / o </> _1 , : <j>(U) — > R". We assume that 



H{p) = H(f(p)) 



(3.1) 
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Fig. 3.1. A map f which satisfies covering conditions. The set is contracted in coordinate y 
and expanded in coordinate x. Note that in the 9 coordinate the set may be simply shifted, expanded 
or contracted, just as long as conditions (3. 3),..., (3. 7) are satisfied. 

for all p £ and also that for some v > condition (2.5) holds. 
We introduce the following sets 

= B? x 8BI x B r s , (3.2) 
D+ = B? x & u x dBl 

We now introduce a number of definitions. The first is a definition of a covering 
relation. 

Definition 3.1. We say that a map f : U — > U satisfies covering conditions in 
Dif 

0, (3.3) 

0, (3.4) 

0, (3.5) 

0, (3.6) 

and for any point p £ B^ x {0}, 

n(x,y)U(j?), H*,v)U l (P) e int x B r s ) • (3.7) 

Conditions (3.5) and (3.4) mean that, in the y (stable) projection, contracts the 
set Dfj, strictly inside B r s . Conditions (3.6) and (3.3) mean that, in the x (unstable) 
projection, expands the set strictly outside B r s . The final assumption (3.7) 
is needed to ensure that the image of D$ by f$ intersects D^. Without assumption 
(3.7), all other assumptions (3.3),. . . ,(3.6) could easily follow from having image of D 
disjoint with D. 

Covering relations are tools which can be used to ensure existence of an invariant 
set in D. To prove that this set is a manifold we shall need additional assumptions. 
These shall be expressed by "cone conditions". To introduce these conditions, first 
we need some notations. 

Let Q h , Q v : R c x R s x R" ->Ebc functions defined by 

Q h (6,x,y) = - 7 ||0|| 2 + a h \\x\\ 2 - (3 h \\y\\ 2 , (3.8) 
Q v (0,x,y) = - 7 \\0\\ 2 - a v \\x\\ 2 + p v \\y\\ 2 , 

with j,ah,ot v , /3fc, f3 v > and 

a h > a v and [} v > f3 h . (3.9) 

Definition 3.2. We say that a map f : U — >• U satisfies cone conditions in D if 
there exists an m > 1 such that 



7r x (U(D^))nB r u = 
n y (U(Dt))n(m s \B r s ) = 

nx(fZHD4>))n(W\K) = 
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Fig. 3.2. An example of a function f, which satisfies cone conditions: A. Two points pi,P2 
for which Qh{pi — P2) = c > 0. B. Difference of the images of the points lie on a cone Qh(ft/>(pi) ~ 
fdipz)) > mc - Similar condition (but with reversed roles of the x and y coordinates) needs to hold 
for the inverse map. 




Fig. 3.3. A horizontal disc in D^. 

1. for any two points pi,P2 €E satisfying p\ 7^ P2 and Qh{pi — P2) > we 
have 

Qh(U(pi) - U(P2)) > mQhipx - p 2 ), (3.10) 

2. for any two points Pi,p 2 € D4, satisfying pi ^ p 2 and Q v {p\ — p 2 ) > we 
have 

Qvif^iPi) - f^(P2)) > mQ v ( Pl -p 2 ). (3.11) 

Definition 3.2 intuitively states that if we have two points that lie horizontally 
with respect to each other, then their images are going to be pulled apart in the 
horizontal, x coordinate (see Figure 3.2). If on the other have we have two points 
that lie vertically with respect to each other, then their pre-images are going to be 
pulled apart in the vertical, y coordinate. 

We now give definitions of horizontal discs and vertical discs. These will be 
building blocks in our construction of invariant manifolds. 

Definition 3.3. We say that a continuous monomorphism h : W u — > is a 
horizontal disc if ir x h(x) — x and for any X\,X2 G B r u 

Q h (h( Xl ) - h{x 2 )) > 0. (3.12) 

Thus, to any point x in the graph h(x) we can attach a horizontal cone, so that the 
graph always remains entirely inside the cone (see Figure 3.3). 

Definition 3.4. We say that a continuous monomorphism v : B 7 S — > is a 
vertical disc if n y v(y) = y and for any y\,yi € B r s 



Q,(v(yi)-v(y 2 )) >0. 
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Thus, to any point y in the graph v(y) we can attach a vertical cone, so that the graph 
always remains entirely inside the cone. 

The following lemma is a key auxiliary result for the proof of Theorem 3.7, which 
is the main result of this section. Roughly speaking, it states that under appropriate 
conditions, an image of a horizontal disc is a horizontal disc. 

Lemma 3.5. Let hi be a horizontal disc. If f satisfies covering and cone condi- 
tions in D, then there exists a horizontal disc h 2 such that 

{p : TT XtV p G K x Bl) n / (hiOO = h 2 (BQ. 

Moreover, if H (</> -1 (hi(i?^))) < h, and for any p € such that H{(f>~ 1 {p)) < h we 
have 



the 



MU(p)) e Bf, (3.13) 



h 2 K)c^ and H{4>- 1 {MB r u )))<h. 



Proof. Without loss of generality we can assume that <f> is equal to identity. Thus 
we can set D<p = D and f^ = f. 

Let h be any horizontal disc, then by (3.8), (3.12) and (3.10) for X\ 7^ x 2 

a h IK/Oi^i)) - ir x f(h(x 2 ))\\ 2 > Q h (f(h( Xl )) - f(h(x 2 ))) (3.14) 

> mQ h (Yi(xi) - h(x 2 )) 
>0, 

which means that n x o f o h is injective. 

Let us define a function F : B r u — >■ E" as follows 

F(x) :=^(/(hx(x))). 

We shall first show that there exists an xo £ B u such that F(xq) g B r u . Using 
notations hi(x) = (Jiq(x), x, h y (x)) we can define a family of horizontal discs h Q (x) = 
(ahe(x) , x, ahy(x)) . We define a function I : [0, 1] X BJ -> W u as 

l(a,x) := 7r x o fo h a (x). 

By (3.2) and (3.3), since h a {dB r u ) C D~ , for any a E [0, 1] we have l(a, dB r u )^B r a = 0. 
Since, as shown at the beginning of the proof, 

l(a,-) := tt x o f o h a : B r u ^ K u 

is a continuous monomorphism, we either have l(a, B^nB^ = or B r u C int(Z(a, B r u )}. 
This also means that 

inf{||Z(o!,0) -x\ : x e dB r u } > 0, 
and thus the function 6 : [0, 1] — > K defined as 



6{a) := 



Z(a,0)e£?; 

1 l(a,0)$BZ, 
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is continuous. We have 

ir x h a = Q (0) = 0, 
7ryh a= o(0) = 0, 

so condition (3.7) implies 1(0, 0) = ir x o /(h a =o(0)) G -B£, hence (5(0) = 0. Suppose, 
to obtain a contradiction, that F(a;) ^ for all x E B^. This would mean that in 
particular F(0) = 1(1,0) £ B r u , hence 6(1) = 1. This contradicts the fact that 5(0) = 
and <5 is continuous. 

We have shown that there exists an xq £ B r u such that F(xo) € B r u . From (3.3) 
follows that F(dB r u ) n = 0. We also know that F = ir x ° / ° hi is continuous and 
injective. Putting these facts together gives B r u C F(.B£). This means that for any 
v G B£ there exists a unique x = x(w) € B r u such that .F(x) = v. We define 

h 2 (u) = (h 2 .e(v),v,h 2 , y (v)) := (tt 6 o f oh. 1 (x(v)),v,7r y o f ohi(x(v))). 

For any v\ ^ v 2 , v\,v 2 & B u , by (3.12) and (3.10) we have 

Q h (h 2 (t>i) - h 2 («a)) = Q h (/ o h x (x(«i)) - / o ht(x(v 2 ))) (3.15) 

> mQ h (h 1 (x(v 1 )) - hi(x(u 2 ))) 

> 0. 

Since Q h (h 2 (ui) - h 2 (u 2 )) > 0, 

"ft IK ~ Wall > Ph \\h 2 , y (v\) ~ h 2 ^ y (v 2 )\\ 2 + 7 \\h 2fi (v x ) - h 2 , e (v 2 )\\ 2 
> mm(/3 h> i) \\(h 2 ,e, h 2 , y )(v 1 ) - (h 2 . e ,h 2 ^)(v 2 )\\ 2 , 

and therefore h 2 is continuous. 

Finally let us note that (3.1) and h 2 (u) = / o h.\(x(v)) implies H(h 2 (B£)) = 
H(hi(5D) < h. This by (3.13) implies that h 2 (5;) C D. □ 

Next lemma follows from mirror arguments. 

Lemma 3.6. Let vi be a vertical disc. If f satisfies covering and cone conditions 
in D, then there exists a vertical disc v 2 such that 

{p ■■ **, v p e K x B r s } n U(MB r s )) = v 2 (B r s ). 

Moreover, if H(4>~ l (vi(B r s ))) < h, and for any p £ such that H(4>~ 1 (p)) < h we 
have 

veU^WeB?, (3.16) 

then 

v 2 (B r s )cD^ and H^ 1 (v 2 (B 7 s ))) < h. 

We are now ready to state our main result for maps, which will be the main tool 
for the proof of Theorem 2.4. 

Theorem 3.7. If f satisfies covering and cone conditions in D, and in addition 
for any p € with H((f>~ 1 (p)) < h we have 

*oU(p)eB* and votftoeB?, (3.17) 
then there exists a C° function x ■ B^~ v ~^ such that 
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1. For any 9 G B^ v we have ngx(O) = 9 and 

/ n (0 _1 (x(0))) G aZZ n G Z. 

2. ///or some p G 4>~ 1 (B^~ V x x B£) we ftawe 

f n (p)eD forallneZ, 

then there exists a 9 £ B^~ v such that p — cj)^ 1 (x(9)). 
Proof. Without loss of generality we assume that <f> is equal to identity, which 
means that = D and f$ = f- 

Let #o G B^~ v and y$ £ B r s . Let hi : B r u ->Dbea horizontal disc defined by 

hi(x) = (9 ,x,y ). 

Clearly hi satisfies cone conditions and also by (2.5), (</) _1 (hi(i?^))) < h. Applying 
inductively Lemma 3.5 we obtain a sequence of horizontal discs hi, h 2 , . . . such that 

/(h,_ 1 K))nfl = h,(fl;) and H(hi(B r a )) < h. 

This by compactness of B r u ensures existence of a point Xq G B r u such that for all 
n G N 

fOn^GD. (3.18) 

Suppose that we have two points x\ and x\ which satisfy (3.18). Then by (3.15) we 
have 

a h r 2 > a h \\n x (/"(h^)) - /"(h^)))!) 2 
>Q,(/"(h 1 (4))-/"(h 1 (^))) 

> mQ h (r^^ixl)) - (3-19) 

>m"Q h (hi^-hi^)), 

which since m > 1 cannot hold for all n. This means that functions W cs : B^~ v x 
->• D, w cs : 3*-° x B£ given as 

W cs (0 o ,2/o) = (0o,w cs (9 o ,y Q ),y Q ) := (0o,4,yo), 

are properly defined. Note that by a similar argument to (3.19), for any (#i,yi) 7^ 
(02, y 2) we must have 

> Q h (W cs (6 2 ,y 2 ) - W cs (9 1 ,y 1 )). (3.20) 

This gives 

max( 7 ,/3 ft ) ||(0 2 ,y 2 ) - (0i,2/i)|| 2 > 7 ||0 a - 9^\ 2 + (3 h \\y 2 - yi \\ 2 

>a h \\w cs (6 2 ,y 2 )-w cs (6 2 ,y 2 )\\\ 

which means that w cs is Lipschitz with a constant 

L 8 = J max(7 ' ft) . (3.21) 
V a/j 
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Mirror arguments, involving Lemma 3.6, give existence of functions W cu : B^ v x 
B r u -> D, w cu : B?-° x B r u -> B r s 

W cu (9 > x) = (6,x,w cu (9,y)), 

such that for any point (9, x) £ B^~ v x B r u and all n e N 

f- n (W cu (6,x)) e D. 

Also w cu is Lipschitz with a constant 



Lu = ^^f^- (3.22) 

We shall show that for any 9 £ B? ~ v the sets W CS {9,B 7 S ) and Vy CM (6»,P£) inter- 
sect. Let us define Pe : B r u x P^ — » P£ x PJ as 

Pe^y) := (7T a W CS (0,y), % T^(0,x)). 

Since Pg is continuous, by the Brouwer fixed point theorem there exists an (#0,2/0) 
such that Pg(x ,y ) — (xo,yo) . This means that 

W cs (9,y ) = (9,w cs (8,y ),yo) = (9,x ,w cu (9,y )) = W cu (9,x ). 

Now we shall show that for any given 9 £ Bf~ v there exists only a single point 
of such intersection. Suppose that for some 9 G B^~ v there exist (xi,yi) , ( 0:2,3/2) G 
P£ x P£, (xi,yi) 7^ (x 2 ,y 2 ) such that 

VF cs (0, yi ) = W cu (9, Xl ) and VF cs (0, y 2 ) = W cu (9, x 2 ). 

We would then have W cs (9, y rn ) = W cu (9, x m ) = (9, x m , y m ) for m = 1, 2. 
From (3.20) follows that 

> Q h (W cs (9, Vl ) - W cs (9, y 2 )) = Q h {{9, Xl , Vl ) - (9, x 2 , y 2 )) , 
and by mirror argument 

> Q v (W cu ( Xl ,X) - W cu (x 2 ,X)) = Q v {(9, Xl , yi ) - (9,x 2 ,y 2 )) 

which implies that 

> (a h - a v ) \\xi - x 2 \\ 2 + (f3 v - f3 h ) \\yi - y 2 \\ 2 , 

which contradicts (3.9). 

We now define x{9) = (9,Xx, v (9)) ■= (9,x ,y ) for x = x (9), y = yo(9) such 
that W cs (9,yo) — W cu (9,Xo)- By previous arguments we know that x is properly 
defined. We need to show continuity. Let us take any 9i,9 2 € B^~ v - From (3.20) 
follows that 

Qh (x(9i) - = Qh (W cs (9 1 ,y (9 1 )) - W cs (9 2 ,y (9 2 ))) < 0, (3.23) 

and by mirror argument 

Qv (x(9i) ~ X{9 2 )) = Q v (W cu (9 1 ,x (9 1 )) - W cu (9 2 ,x (9 2 ))) < 0. (3.24) 
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From (3.23), (3.24) follows that 

an \\xo (#i) - x (6 2 ) 
-a v \\x Q (00 - x a (6 2 )\\ z + (i v \\y (0^ - y (0 2 )\\' < 7 - 6 2 \ 



a h \\x (0x) - x (9 2 )\\ 2 - (3 h \\y (0i) - y (0 2 )\\ 2 <j\\0 1 - 6 2 f 



K - a v ) \\x Q (00 - x (0 2 )|| 2 + (& - h ) \\ yo (00 - y (0 2 )|| 2 < 2 7 \\6 1 - 2 || 2 . 
which gives 



||Xx,„ ~ X*,v (fc)ll < 4/ -r-7 ^ ttt -O2W, (3.25) 

and by (3.9) implies Lipschitz bounds for \x,y and continuity of \- D 

4. Proof of the main theorem. In this section we shall show that assumptions 
of Theorem 2.4 imply that a map induced as a shift along a trajectory of the flow of 
(2.1) for sufficiently small time satisfies covering and cone conditions. This will allow 
us to apply Theorem 3.7 to prove Theorem 2.4. 

We start with a lemma which shows that assumptions of Theorem 2.4 imply 
covering conditions for a shift along the trajectory of (2.1). 

Lemma 4.1. Assume that for any p 6 BfJ~ v x {0} x {0} assumptions (2.20), 
(2.21), (2.19) of Theorem 2.4 hold, then for sufficiently small r > and all t <G (0,t] 
a function 

/(x) :=*(t,x) 

satisfies covering conditions. 

Proof. Without loss of generality we assume that (j> = id. Let q € x {0} x 
{0} n N p . By (2.19), for sufficiently small t 

\\n x $(t, q)\\ = n x f 3(0, g) + ^U(0, q)t + o(t) 
= ||0 + t7T^(g)+o(i)|| 

<\t\rE u . (4.1) 
Analogous computation yields 

\\n v *(t,q)\\<\t\rE.. (4.2) 
In later parts of the proof we shall use the fact that for any q\,q 2 € N„ 



F{qi) - F(q 2 ) = [ dF (g 2 + s (gi - g 2 )) rfs ( 9l - g 2 ) • 
Jo 

Now we shall prove (3.3). Let q = (0, a;, y) € D7 n A^, which means that 



(4.3) 
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Using |$(*,g)| t=0 = F(q), $(0,g) = q, and (4.3) we have 
| IK, ($(i,g)-$ ( ^(^o,0)))|| 2 | t= o 

- j t (ir x g) - *(t, (0, 0, 0))) T tt, (*(t, g) - $(*, (9, 0, 0)))) ^ 
= 2ir x (q - (9, 0, 0)) T n x (F(q) - F (9, 0, 0)) 
= 2x t -k x (J dF (9,sx,sy)ds(0,x,y)^j 
= 2x T (Ax + e u y) , 

where 

A = j ~ix~ ^' SXl Sy ^ dSl 6u ^J ~^y^~ ^' SX ' ^ dS ' 
From bounds (2.9) and (2.12) we thus obtain 

j t \\tt x (<£(*, g) - S(t, (0, 0, 0)))|| 2 | t=0 > 2 (rV - ||x|| ||e u || ||y||) > 2r 2 (S u - e u ) . 

(4.4) 

Using the same arguments we can also show that for any q = (8, x, y) £ N p 

| K mt,q) - $(*, (0,O,O)))|| 2 | t=0 < 2 ||y|| (e s r - ||y|| 5 s ) (4.5) 
Combining (4.1) (4.4) and (2.20), for sufficiently small t > gives 

11^/(^11=117^(^)11 

> \\tt x ($(*, g) - $(t, (0, 0, 0)))|| - |kx* (*, (0, 0, 0))|| 

> \fh* ($(t,g)-$(i,(0,O,O)))|| 2 - ir£ u (4.6) 

> ^ |K (*(0, q) - *(0, (0, 0, 0)))|| 2 + <2 (<5" - e u r 2 ) - ir£ u 
= Vr 2 + t2r 2 (<S« - e u ) - trE u 

> r. 

This establishes (3.3). Now we shall show (3.5). For any g = (9,x,y) <G and 
sufficiently small t > 0, analogous derivation to (4.6) (for these computations we use 
estimates (4.2), (4.5)) give 

K/(?)ll < \/\\v\\ 2 +t2\\y\\(e a r-\\y\\5») + trE s . (4.7) 

Since \\y\\ < r by (2.21), for sufficiently small t > 0, inequality (4.7) implies that 
ll 7r y/( ( i , )ll < r an d hence establishes (3.5). 

Proof of (3.4) and (3.6) follows from analogous arguments with t < 0. 

Conditions (3.7) hold for sufficiently small t. This follows from continuity of $(p, t) 
with respect to t since 

/(p) = $(t,p) r i ( J3 ) = $(-i, P ), 
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and for p e B* x {0} x {0} 

n {x<y) ${0,p) = n (x<y) p = (0, 0) G int (B r u x B r s ) . 

□ 

Now we shall show that assumptions of Theorem 2.4 imply cone conditions for a 
shift along trajectory of (2.1). Let us start with a simple technical lemma. 

Lemma 4.2. Let C — {Cij) i . =1 2 3 be a (c + u + s) x (c + u + s) matrix. Assume 
that for a,-, bi g R, i = 1, 2, 3 we /iawe 

inf{sf : H^ll = 1} > a;, /or* = 1,2,3, (4.8) 
sup{a;f Cii^i : INI = 1} < &< /or » = 1, 2,3, (4.9) 

</ien for any x = (x\, X2, £3) € u c +'"+ s 

x T Cx > {at - ci) \\ Xl \\ 2 + (a 2 - c 2 ) ||a; 2 || 2 + (03 - c 3 ) ||x 3 || 2 , (4.10) 
o^Cz < (&! + ci) Hxill 2 + (62 + c 2 ) ||x 2 || 2 + (63 + c 3 ) ||x 3 || 2 , (4.11) 

where 

ci = i(||c 2 i|| + ||a3i|| + ||c 12 || + ||c 13 ||), 

C2 = ^(||C 2 i|| + ||C 2 3|| + ||C 12 || + ||C 32 ||), 
C8=^(||C3 1 || + ||C 23 || + ||C 13 || + ||C32||). 

Proof. The estimate (4.10) follows by direct computation from (4.8) and the fact 
that for any i, j 

±2qJC Jl q l >-\\C Jl \\(\\q ] \\ 2 + \\q l \\ 2 ). 
Similarly (4.11) follows from (4.9) and 

±2qfC ji q i <\\C ji \\(\\q j \\ 2 + \\q i \\ 2 ). 

□ 

Let Ik denote a k x k identity matrix. Let 

Qi = diag(-7/ c , a h I u , -Phis), 
Q-2 = diag(-7/ c , -a v I u , f3 v I s ), 

be matrices associated with quadratic forms Qh and Q v respectively. Now we are 
ready to prove that assumptions of Theorem 2.4 imply cone conditions for a time 
shift along a trajectory map. 

Lemma 4.3. Assume that for any p € B^~ v x {0} x {0} assumption (2.17) of 
Theorem 2.4 holds, then for sufficiently small r > and all t £ (0, t] a function 

/(x) :=$(t,x) 

satisfies cone conditions with a coefficient m — 1 + th, with some constant h > 0.. 
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Proof. Let pi,p 2 € be such that Pi = (9i,Xi,yi) for i = 1,2, pi ^ p 2 and 
<2h(Pi -P2) > 0. Let p = (0i,O,O) £Bf x {0} x {0}. Condition (2.14) implies that 
Pi,P2 € iV p . We compute 

^ (($(t,pi) - $(i, P2 )) T Qi(<f (t,pi) - *(t,pa))) k=o 

= 2(p 1 -p 2 ) T Q 1 (F( Pl )-F( P2 )) (4.12) 

= 2(p 1 -p 2 ) T Q 1 B(p 1 -p 2 ), 

where 

B = / dF(p 2 + t(pi - p 2 ))dt e [dF(Np)}. 
Jq 

For C = QxB, from (2.11), (2.9), (2.10) we have 

mi{xfCiiZi : ||xi|| =1} > - 7 c u , 

mi{x%C 22 x 2 : \\x 2 \\ = 1} > a h S u , (4.13) 
inf^C^a : \\x 3 \\ - l}>PhS s . 

Using (4.10) from Lemma 4.2 with (4.13) and (2.12), for 4° rw , n f ° rw , n f s mw given by 
(2.15) and m G (max(4° rw , 4° rw ), < orw ) we have 

a; T Cx > -4°™ 7 ||zi|| 2 + < orw a h ||x 2 || 2 - p h \\x 3 f 
> Ml (-7 ll^ill 2 + a h \\x2W 2 - Ph \\xz\\ 2 ^j 

= nix T Qix. (4-14) 

The constant /Ui g (max(«;[° rw , w ), k^ ™) can be chosen to be greater than zero 
thanks to assumption (2.17). This means that by (4.12) and (4.14) 

- ((${t, Pl ) - $(t,p 2 )fQ 1 (<f>(t,p 1 ) - Ht, P2 ))) \ t=0 > 2fi 2 Q h ( Pl -P2). 
For sufficiently small r > and t <E (0, r) we therefore have 

Qh(f(Pi) - /(pa)) = Qh(*(*,Pi) ~ *(*,!*)) 

= Qh(pi ~ P2) + t—Q h mt,px) - $(t,p 2 ))|t=o + o(t) 
> (l + t2/ii)Q fc (pi-pa), 

which establishes (3.10) with m = 1 + i2/.«i > 1. 

The proof of (3.11) is obtained analogously with m = 1 + i2/x 2 > 1 for some 
M2 < 0, /i 2 € (Kl? ack ,min(^ ack ,Kjf ck )) , with negative time t < 0. 

So far the entire argument was done for points in AT p . We can choose h p = 
min{2 |/Zi| ,2 |/i 2 |} so that (3.10) and (3.11) hold for any pi,p 2 G iV„ with a constant 
m = 1 + |t| ftp. By compactness of we can now choose a h > such that (3.10) 
and (3.11) hold with a constant to = 1 + |i| /i for all pi,p2 € D 

We are now ready for the proof of our main result. 
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Proof. [Proof of Theorem 2.4]By Lemmas 4.1 and 4.3 we know that assumptions 
of Theorem 2.4 imply cone and covering conditions for a map induced by the flow by 
a small time shift. Now we just need to show that for a map 

/(x) :=$(*, x), 

with sufficiently small t > 0, for any p £ with H((f)^ 1 (p)) < h we have (3.17). This 
follows from (2.5) and continuity of $(t,x) with respect to t. The claim now follows 
from Theorem 3.7. □ 

By applying Theorem 3.7 in our proof of Theorem 2.4 we have established more 
than just continuity of our center manifold. We have also obtained existence of its 
stable and unstable manifolds, together with explicit Lipschitz type bounds on their 
slopes. This is summarised in the following corollary. 

Corollary 4.4. During the course of the proof of Theorem 3.7 we have shown 
that in local coordinates given by </> the stable, unstable and center manifolds obtained 
by our argument are given in terms of functions 

W cs . gR-v x S r D ^ 
W cu . gR-v x £r ^ D ^ 

X : B?~ v -> D+, 

respectively. We have also shown that these functions are of the form 

W cs (9,y) = (e,w cs (6,y),y), 
W cu (e,x) = (9,x,w™(9 ,y)), 

x(0) = (0,x*,y(0)), 

with functions w cs : B^' v x B r s -> B r u , w cu : B^ v x B r u -> B r s and Xx, y ■ B^ v -> 
B 7 U x B r s by (3.21), (3.22) and (3.25) satisfying Lipschitz conditions with constants 

/ max(7,/?h) 
/ max(7, a v ) 

V Jv ' 

/ 2 7 
y mm{a h - a v ,0 v - p h ) ' 

Thus our method gives explicit Lipschitz type bounds for our invariant manifolds of 
(2-1)- ' 

5. Centre manifold around L\ in the Restricted Three body problem. 

In the following we specialise our study to the center manifold of the equilibrium point 
L\ in the restricted three body problem, or RTBP for short. 

Section 5.1 describes the RTBP and presents its equations of motion and specifies 
the equilibrium point L\ around which we shall later prove existence of the center 
manifold. A general reference for this section is Szebehely's book [S]. Section 5.2 
constructs "aligned coordinates" (described in Section 2.1) around L\ in the RTBP 
using a suitable normal form procedure. A general reference for this section is the 



L s = 
L u = 
L c = 
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Fig. 5.1. Notation for the rotating system of coordinates with origin at the center of mass. 
The Sun has the mass 1 — fj. and is fixed at Pi = (fJ.,0). The Earth has the mass fi is fixed at 
P2 = (fi — 1,0). The third massless particle moves in the XY plane. 



paper by Jorba [ ] on computation of normal forms with application to the RTBP. In 
Section 5.3 we show how normal forms can be used to obtain a very accurate numerical 
estimate on where the centre manifold is positioned. In Section 5.4 we apply Theorem 
2.4 to obtain a rigorous enclosure of the centre manifold. 

5.1. Restricted Three Body Problem. The problem is defined as follows: 
two main bodies rotate in the plane about their common center of mass on circular 
orbits under their mutual gravitational influence. A third body moves in the same 
plane of motion as the two main bodies, attracted by the gravitation of previous two 
but not influencing their motion. The problem is to describe the motion of the third 
body. 

Usually, the two rotating bodies are called the primaries. We will consider as 
primaries the Sun and the Earth. The third body can be regarded as a satellite or a 
spaceship of negligible mass. 

We use a rotating system of coordinates centred at the center of mass. The plane 
X, Y rotates with the primaries. The primaries are on the X axis, the Y axis is 
perpendicular to the X axis and contained in the plane of rotation. 

We rescale the masses fti and pi2 of the primaries so that they satisfy the relation 
Ml + M2 — 1- After such rescaling the distance between the primaries is 1. (See 
Szebehelly [S], section 1.5). 

Let the smaller mass be pt 2 = pi = 3.040423398444176 x 10~ 6 and the larger one 
be Hi — 1 — /i, corresponding to the values of the Earth and the Sun respectively. We 
use a convention in which in the rotating coordinates the Sun is located to the right 
of the origin at Pi = (fi, 0), and the Earth is located to the left at P 2 = (fi — 1, 0). 

The equations of motion of the third body are 

x -2Y = n x , 

Y + 2X = Q Y , 



where 



n= l -{x 2 + Y 2 )+ l -^ + ^. 

2 n r 2 
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and r\ , ri denote the distances from the third body to the larger and the smaller 
primary, respectively (see Figure 5.1) 

r 2 = (X-^) 2 + Y 2 , 
r\ = (X-fi + lf + Y 2 . 

These equations have an integral of motion [S] called the Jacobi integral 

C = 20- (X 2 + Y 2 ). 

The equations of motion take Hamiltonian form if we consider positions X, Y 
and momenta Px = X — Y, Py = Y + X. The Hamiltonian is 

H= 1 -{P 2 X +P 2 )+YP X -XP Y - 1 -^-^, (5.2) 
2 n r 2 

with the vector field given by 

F = JVH, 



J = 



( id \ 








I -id 


, id = 


1 


i 



The Hamiltonian and the Jacobi integral are simply related by H = — y . 

Due to the Hamiltonian integral, the dimensionality of the space can be re- 
duced by one. Trajectories of equations (5.1) stay on the energy surface M given 
by H(X,Y, Px , Py) — h — constant, a 3-dimensional submanifold of M 4 . Equiva- 
lently, M is the level surface 

M = {C(X,Y,X,Y) = c=-2h} (5.3) 

of the Jacobi integral. 

The restricted three body problem in a rotating frame, described by equations (5.1), 
has five equilibrium points (see [S]). Three of them, denoted Li,L^ and L3, lie on 
the X axis and are usually called the 'collinear' equilibrium points (see Figure 5.1). 
Notice that we denote L\ the interior collinear point, located between the primaries. 

At this point we would like to make it clear that in this paper we focus only on the 
equilibrium point L\, though other collinear equilibria points could be investigated in 
the same manner. 

The Jacobian of the vector field at L\ has two real and two purely imaginary 
eigenvalues. Since the three body problem is Hamiltonian in can be shown by the 
Lyapunov-Moser theorem [M] that in a sufficiently small neighbourhood of L\ there 
exists a family of periodic orbits which is parameterised by energy. This family of 
orbits forms a center manifold. Our aim shall be to prove the existence of this manifold 
in a neighbourhood which is far from L\, As mentioned before, close to L\ the existence 
of this manifold follows from the center manifold theorem (or in this case also from 
the Lyapunov-Moser theorem). The hard task is to prove its existence far from the 
equilibrium point. 

Remark 5.1. Since the center manifold around L\ is foliated by periodic orbits 
it has to be identical to the invariant manifold obtained through Theorem 2.4 due 
to point 2. of the theorem. The Lyapunov-Moser theorem ensures the existence of 
periodic orbits locally. In such local domain we are guaranteed that the manifold 
X from Theorem 2.4 inherits all regularity properties which follow from the center 
manifold theorem. Outside of this domain Theorem 2.4 establishes only Lipschitz 
continuity of X- 
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5.2. Normal Form. The linearised dynamics around the equilibrium point L\ 
is of type saddle x centre for all values of /i. In this section we use a normal form 
procedure to approximate the nonlinear dynamics locally around L\. 

For the purpose of this paper, the normal form coordinates will be used precisely 
as the well-aligned coordinates described in section 2.1. 

The goal of the normal form procedure is to simplify the Taylor expansion of the 
Hamiltonian around the equilibrium point using canonical, near-identity changes of 
variables. This procedure is carried up to a given (finite) degree in the expansion. 
The resulting Hamiltonian is then truncated to (finite) degree. Such Hamiltonian is 
said to be in normal form. 

We compute a normal form expansion that is as simple as possible, i.e. one that 
has the minimum number of monomials. This is sometimes called a full, or complete, 
normal form. The equations of motion corresponding to the truncated normal form 
can be integrated exactly. As a result, locally the normal form gives a very accurate 
approximation of the dynamics. 

In particular, here we use the normal form to approximate the local center man- 
ifold by a 1-parameter family of periodic orbits with increasing energy. 

The normal form construction proceeds in three steps. First we perform some con- 
venient translation and scaling of coordinates, and expand the Hamiltonian around L\ 
as a power series. Then we make a linear change of coordinates to put the quadratic 
part of the Hamiltonian in a simple form, which diagonalises the linear part of equa- 
tions of motion. Finally we use the so-called Lie series method to perform a sequence 
of canonical, near-identity transformations that simplify nonlinear terms in the Hamil- 
tonian of successively higher degree. 

The transformation to well-aligned coordinates (j> : U — > 4>{U) C W L is the com- 
position of all the transformations performed during these three steps. 

A similar full normal form expansion has been used for the spatial RTBP in a 
previous paper [DMRj. We refer to the previous paper for the fine details of the 
normal form construction, which will be left out of the current paper. 

5.2.1. Hamiltonian expansion. We start by writing the Hamiltonian (5.2) as 
a power scries expansion around the equilibrium point L±. First we translate the origin 
of coordinates to the equilibrium point. In order to have good numerical properties for 
the Taylor coefficients, it is also convenient to scale coordinates [ ]. The translation 
and scaling are given by 



where 7 is the distance from L\ to its closest primary (the Earth). 

Since scalings are not canonical transformations, we apply this change of coordi- 
nates to the equations of motion, to obtain 



X = -jx + /1 - 1 + 7, Y = -7?/, 



(5.4) 



x — 2y = 
y + 2x = Q. 



X 



where 



H - 1 + 7 



x + 



1_ / 1 - n ^ 



) 



7 



7 3 V r i r 2 



and ri , T2 denote the (scaled) distances from the third body to the larger and the 
smaller primary, respectively. 
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Defining p x — x — y, p y — y + x, the libration-point centred equations of mo- 
tion (5.5) are Hamiltonian, with Hamiltonian function 



# = + Py> + fP* ~ xp y + 



1 / 1 - n 



7 7 \ r i 

Our first change of coordinates can therefore be summarised as R : R 4 — > M.' 

(X,Y 1 P x ,P Y )=R(x,y lPx ,p y ) 

= (-jx + n - 1 + 7, -72/, -7p x , -7P„ + M - 1 + 7) 

The Hamiltonian is then rewritten in the form [J, JM] 

1 



H = ^(Px+Py) + VPx ~ Xp y 



(5.6) 



(5.7) 



J2^)P n Pn(^, (5- 



where P n is the n-th Legendre polynomial, and the coefficients c„(/i) are given by 



1 ( , sr, (1 - M)7 n+1 
c„( A1 ) = - /i+(-i)" 1 



7 

This expansion holds when p < min(|Pi|, IP2I 
centred at L\ that extends up to the Earth. 



(1 - 7)"+ 1 
= IP2I = 1 5 i-e. it is valid in a ball 



5.2.2. Linear changes of coordinates. Now we transform the linear part of 
the system into Jordan form, which is convenient for the normal form procedure. This 
particular transformation is derived in [J, JM], for instance. 

Consider the quadratic part P2 of the Hamiltonian (5.8), 



H 2 = \{pi+ Py) + VPx ~ xp y - c 2 x z + ^y z 



(5.9) 



which corresponds to the linearisation of the equations of motion, ft is well-known [JM] 
that the linearised system has eigenvalues of the form ±A, where A, v are real 
and positive. 

One can find ([JM] section 2.f) a symplectic linear change of variables 



2 A 



( 

X 2 -2c 2 -l 



c 



-2A 

2 Sl 
A 2 -2c 2 -l 

si si 

A 2 +2c 2 + l A 2 +2c 2 + l 

si si 

A 3 -t-(l-2c 2 )A ~A 3 -(l-2e 2 )A 

si si 





-v 2 -2c 2 -l 

2 32 
S2 





2r 
S2 





-V i + (l-2c 2 )v 

S2 



j 



where 



1 = ^2A ((4 + 3c 2 ) A 2 + 4 + 5c 2 — 6cf), 

2 = ^((4 + 3c 2 )« 2 -4-5c2+6c|), 



that puts the linear terms of the vector field at L\ into a Jordan form. This means 
that the change from position-momenta to new variables {x\, y\, X2, y%) € R 4 , 



{x,y,Px,p y ) = C(x 1 ,y 1 ,x 2 ,y 2 ), 



(5.10) 
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casts the quadratic part of the Hamiltonian into 

H 2 = Xx iyi + ^(xl + y 2 2 ). (5.11) 

The linear equations of motion (£, y) — A(x, y) associated to (5.11) decouple into 
a hyperbolic and a center part 

{ii,vi)=A h (xi,yi) (5.12a) 
(x 2 ,y 2 ) = A c {x 2 ,y 2 ), (5.12b) 

with 

A h = (J _° A ) and A c = ( ^ V Q 

Notice that the matrix A of the linear equations (5.12) is in block-diagonal form. 
It is convenient to diagonalise the matrix A over C. Consider the symplectic change 
T~ 4 : M 4 — > C 4 to complex variables (q,p) = (qi,Pi, qi,Pi) € C 4 

(<?1jP1j92,P2) = T~ 1 (x 1 ,y 1 ,x 2 ,y 2 ) (5.13) 
^( x 2 — is/a), ^(~ ix2 + y2 ) 
This change casts the quadratic part of the Hamiltonian into 

H 2 = XqiPi + ivq 2 p 2 . (5-14) 
Equivalently, this change carries A to diagonal form: 

T~ l AT = A = diag(A, -A, iv, -iv). 

5.2.3. Nonlinear normal form. Assume that the symplectic linear changes of 
variables (5.10) and (5.13) have been performed in the Hamiltonian expansion (5.8), 
so that the quadratic part H 2 is already in the form (5.14). 

Let us thus write the Hamiltonian as 



H(q,p) = H 2 (q,p) + H 3 (q,p) + H 4 (q,p) + ■■■ (5.15) 

with Hj(q,p) as homogeneous polynomials of degree j in the variables (q,p) G C 4 . 

As shown in a previous paper [DMR], we can remove most monomials in the 
series (5.15) by means of formal coordinate transformations, in order to obtain an 
integrable approximation to the dynamics close to the equilibrium point. 

Proposition 5.2 (Complete normal form around a saddle x centre) . 

[DMR] For any integer N > 3, there exists a neighbourhood IA^ N ' of the origin 
and a near-identity canonical transformation 

T {N) : C 4 D»M ^C 4 (5.16) 

that puts the system (5.15) in normal form up to order N , namely 



:= H o T [N) =H 2 + i?W + K.W 
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where Z^ N ^ is a polynomial of degree N that Poisson- commutes with H 2 

{ZW ff 2 } = 0, 

and is small 

i^ (Ar) (z)i<c JV ||zir +i v«6«w. 

If the elliptic frequencies v, uj are nonresonant to degree N , 

Ci v + c 2 uj^0 V(ci,c 2 )eZ 2 , < |ci +c 2 | < N, 

then in the new coordinates, the truncated Hamiltonian H 2 + Z^ N > depends only on 
the basic invariants 

h = QiPi = XiVi (5.17a) 
x^ -(- y^ 

h = mP2 = <72<Z2 = 2 2 2 ■ (5.17b) 

The equations of motion associated to the truncated normal form H 2 + Z^ N ^ can 
be integrated exactly. 

Remark 5.3. The reminder 1ZS N ' is very small in a small neighbourhood of the 
origin. Hence, close to the origin, the exact solution of the truncated normal form is 
a very accurate approximate solution of the original system H . 

Remark 5.4. Let <f>i,(f>2 be the symplectic conjugate variables to I\,Ii, respec- 
tively. The basic invariant I2 is usually called action variable, and its conjugate vari- 
able cf>2 is usually called angle variable. They are given in polar variables (5.17b). 

We can now write our function <j) for our change into the well aligned coordinates 
(2.3). To do so we compose the inverse transformations given in (5.7), (5.10), (5.13) 
and (5.16) which gives us 

0= (t^Y 1 oT- x oC- l oET x . (5.18) 

Remark 5.5. The above described method of obtaining normal form coordinates 
is performed by passing through complex variables. It is possible though to arrange 
the changes so that the combined change of coordinates (5.18) passes from real to real 
coordinates. The change of coordinates <j) is a high order polynomial. It is possible 
to arrange the normal form change of coordinates so that the coefficient of <f> are real 
(see [ ]). In setting up our change of coordinates for the application of Theorem 2.4 
to the RTBP in Section 5.4 we have adopted such a procedure. 

Remark 5.6. In practice, one usually computes a normal form of degree N = 
16. In our application to the restricted three body problem in Section 5.4 we use a 
normal form of degree N = 4. This turns out to be sufficient, since we investigate a 
relatively close neighbourhood of the invariant point, where degree of order four gives 
us a sufficiently good approximation. 

5.3. Approximating the center manifold in normal form coordinates. 

In normal form coordinates given by (5.18) the Hamiltonian, by Proposition 5.2, is of 
the form 



n (N) =h 2 + zW +n m , {z^ n \h 2 } = o. 



(5.19) 
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In this section we shall show that when we neglect the reminder term , and thus 
consider an approximation of the system, the normal form coordinates given by (5.18) 
give us a very good understanding of where the center manifold is positioned and of 
the dynamics on it. 

Let U be some small neighbourhood of the fixed point (in our discussion for 
the R3BP this will be L\) and let (j>: U — > 4>(U) C K 4 be the transformation to 
normal form coordinates (5.18). Consider the normal form (5.19) up to order N with 
associated equations of motion 

p = F+(p) := JVH (N) (p). (5.20) 

Consider now the truncated normal form up to order N 

HW = H 2 + ZW, 

with associated equations of motion 

p = F+(p) := JVH iN) (p). (5.21) 

Recall that the corresponding linearisation around the origin is (5.12) 

P = Ap, p=(x 1 ,y 1 ,x 2 ,y 2 )eM 4 , (5.22) 

where Xi,yi are the hyperbolic normal form coordinates (5.17a), and #2,2/2 are the 
center normal form coordinates (5.17b). In order to match the notation from section 2, 
let us denote the center normal form coordinates x 2 , y 2 as 9i,9 2 , the unstable normal 
form coordinate X\ as x, and the stable normal form coordinate j/i as y. Note that to 
match the notations we need to swap the order in which the coordinates are written 
out passing from (x\ , y% , x 2 , y 2 ) to (9i,9 2 ,x,y). 

The truncated system has several invariant subspaces. Specifically, the next 
proposition follows from [Mu], section 5.1. 

Proposition 5.7. Let 

E c = {(01,02,0,0): (0 1; 2 ) £M 2 }, (5.23) 
E u = {(0,0,x,0): x e M}, (5.24) 
£ s = {(0,0,0,y): yei}. (5.25) 

Then, E c , E u and E s are invariant subspaces of the flow of . 

Remark 5.8. These subspaces are invariant of the nonlinear truncated system 
(5.21), not just of the linearised system (5.22). It is important to stress here though 
that these subspaces need not be invariant under the full system (5.20). 

Next we claim that E c is approximately equal to the center manifold W c of the 
full system F&. This is formulated in the next proposition which follows from [Mu], 
Section 5.2. 

Proposition 5.9. For each integer r with N < r < 00, there exists a (not 
necessarily unique ) local invariant center manifold W c of F^ of class C r such that 
• W° is expressible as a graph over E c , i.e. there exists a neighbourhood V C E c 
and a map x- V — >• E u ffi E s such that 



W c = {{9 1 ,9 2 ,x,y)£E c ®E u (BE s : {9 1: 9 2 )eV, {x,y) = X {9i,9 2 )}. 
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• W c has N-th order contact with E c , i.e. x an d ^ s derivatives up to order N 
vanish at the origin. 

Hence, in normal form coordinates, the center manifold W c of is approximated 
very accurately (to order N) around the origin by the subspace E°. 

Remark 5.10. When applying Proposition 5.9 we are faced with a problem that 
it is usually very hard to obtain a rigorous bound on the size of the set V . Moreover, 
even though we know that up to order N the derivatives of x vanish at zero, it is 
usually very hard to obtain a rigorous bound on the size of the higher order terms of 
X on the set V, and thus obtain a rigorous bound on the position of the true centre 
manifold. 

Let us now briefly discuss the dynamics of the system (5.21) on E c . To do so we 
shall use the center normal form coordinates (5.17b) in action-angle form, i.e. from 
now on we will use (I, cp) 6ixT for the center part. Proposition 5.2 states that the 
truncated Hamiltonian depends only on the action /, and not on the angle tp. 

Thus the restriction of to its invariant subspace E c is 

1 = 0, ^= JL—=: UJ {I). (5.26) 

The solutions inside E c with initial condition 1(0) = Iq and <p(0) = ifio are I(t) = Iq, 
<p(t) = u}(I )t + <fQ. In the case of the restricted three body problem E c is two 
dimensional and so the dynamics of the truncated system on E c is foliated by invariant 
circles of increasing action /. Notice from equation (5.11) that H grows linearly with 
respect to I (close to the origin), so the invariant circles also have increasing energy 
H. 

The properties discussed above motivate the use of the normal form coordinates 
01, 2 , x, y as the well-aligned coordinates in the sense of section 2.1. They provide a 
good guess on the location of the center manifold (locally around the origin) . Taking 
Bf = {(01,02) e K 2 : ||0i,0 2 || < R} the guess is given by cj>- x {B? x {0}) C M 4 . 
Notice also that the center coordinate I is well-aligned with the energy (in sense of 
(2.5)). Let 

C fl ={(0i,0 2 )€M 2 : ||0i,0 2 || =R} = j(J, <p) € R X T : I = 

be the invariant circle of radius R for the system (5.21). By (5.26) we have i-LS N '(C Rl x 
{0}) < n^{C R2 x {0}) whenever J2 X < R 2 . Hence, given an energy h, we can find 
i?i , i?2 > such that 

u( N \c Ri x {o}) < h < n {N) {c R2 x {o}). 

Taking R\ , i? 2 sufficiently far (in practice they are still close) from one another and 
taking sufficiently small r > 0, since and n w are close, we expect also that 

H {N) {Bf* xB 7 u x B r s ) <h< n (N) {C R2 xB r u x B r s ). 

since U (N) = H o (j)' 1 , this will mean that that the bound (2.5) shall be satisfied. 

5.4. Application of the main theorem to the center manifold around 

L\. In this section we shall show how to apply Theorem 2.4 in practice. 

As described in Section 5.2 the change of coordinates to well aligned coordinates 
can be done using a change to normal coordinates (5.18). We obtain the function 
4> using the algorithm of Jorba [ ]. The algorithm allows us to obtain as a real 
polynomial, passing from K 4 to M. 4 . 
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5.4.1. Methodology. To apply Theorem 2.4 it is enough to derive a rigorous 
bound on the derivative of F^. Let us now outline how such a bound can be obtained. 
Using (2.6), for any p £ I 4 we have 

D (F*(p)) = C 2 0(0- 1 (p))^(^ 1 )b)^(r 1 (p)) (5.27) 

In our computer assisted proof we apply the above formula using an interval-arithmetic- 
based software called CAPD (Computer Assisted Proofs in Dynamics 1 ). This soft- 
ware in particular allows for rigorous-interval-enclosure-based computation of high 
order derivatives of functions on sets. In our application wc obtain a global bound 
for the derivative (2.8) on the entire set D^. To compute [DF^{D^)} applying (5.27) 
we only requite to compute images of functions, derivatives of functions and a second 
derivative on a set D^. All such computations can be performed in CAPD. 

Before specifying the size of the set and giving rigorous-interval-based nu- 
merical results, we have to stress one technical problem encountered when applying 
formula (5.27). We take our change to well aligned coordinates to be a high or- 
der polynomial obtained from non-rigorous computations. To apply formula (5.27) 
directly we would need to know its inverse c\T l . Let us stress that one can not use 
a numerical approximation of an inverse change and use it as 0" 1 (such numerical 
approximate inverse is readily available from algorithms of [J]). To apply (5.27) di- 
rectly one would have to use a rigorous, analytic inverse. Since <fi is a polynomial in 
high dimension and of high order, its analytic inverse is next to impossible to obtain 
in practice. To remedy this problem we slightly modify (5.27). Using the fact that 
D(cf)~ 1 )(p) = (l?0((/!> _1 (p))) we can rewrite (5.27) as 

DF*(p) = J D 2 0(0- 1 (p))(^(0- 1 ( P )))- 1 ^(0- 1 (p)) 

+ D^ip^DF^ip)) (D^-^P)))' 1 ■ 

This in interval arithmetic notation gives us the following formula for the interval 
enclosure of D F^ on some set I C 

[DF+(T)] C ^ 2 0([r 1 (I)])W([r 1 (I)]))- 1 ^([r 1 (I)]) (5.28) 
+ D</>( [4>-\l)] )DF( [cb- 1 (I)] ) (D4>{ [r 1 (I)] 

To compute the right hand side of the above equation there is no need to invert the 
function <f>. It is enough to find a set [(/> -1 (I)] which contains the pre-image of I, i.e. 

r 1 ® c [^(i)] , 

and for this we do not need to compute the inverse function. For a set Bel 4 the 
following lemma can be used to verify that </> -1 (I) C B. 

Lemma 5.11. Let <p :W n E™ be a homeomorphism and let I, B C K™ be two 
sets homeomorphic to n- dimensional balls. If (f>(dB)nl — and for some point p G B 
we have 4>{p) € I then 

Proof. This follows from elementary topological arguments. □ 



http: / /capd. ii.uj.edu. pi 
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To apply the lemma in practice it is convenient to first have a non-rigorous guess 
on the inverse function, let us denote it by <j)~ l . This means that 

A function cj)^ 1 is readily available from algorithms of Jorba [J]. We can then choose 
A > 1 and set B = A[0 _1 (I)] (in our application we choose A = 3 which we find is 
large enough for our problem). Then we divide the boundary OB into smaller sets and 
verify that the image by of each smaller set is disconnected with I. We also check 
that for the middle point p in B we have <f>(p) £ I. This by Lemma 5.11 guarantees 
that _1 (I) C B. 

Remark 5.12. Once a set B such that _1 (I) C B is found, there is a useful 
trick that can be used to refine this initial guess on the pre-image. One can take a very 
small set Iq C I and using Lemma 5.11 find a small set Bq such that 4>~ 1 (lo) C Bq. 
The set [</> -1 (I)] can then be chosen as 

[d>-Hl)] =B +[(Dcj)(B))- 1 ] [I-I ]. 

Such choice guarantees that <^~ 1 (I) C [</> -1 (I)] . It is also usually tighter than the 
initial guess B; which is true especially when the function (j) is close to identity. 
Proof. This follows from the mean value theorem. □ 

In a similar fashion to the method from Remark 5.12, to compute the energy for 
a set I C we take some small set Io C I and compute 

wr 1 ®)] c H([r\io)}) + [DHarHm] [ir 1 ®] irHh)]] . (5.29) 

Remark 5.13. When applying above tools to compute [DF^(T)] using (5.28), it 
pays off to use the fact that (f> is composed of linear changes of coordinates, together 
with a nonlinear change T^ N > which is close to identity. Keeping track of both linear 
and nonlinear change allows to tighten the interval bounds of computations. 

To prove the existence of a fixed point (in case of the RTBP we take the point 
Li) inside of our set we use the interval Newton method. 

Theorem 5.14. / ]Let F : R n -> E" be a C 1 function. Let I = njLja^&i], 
Oi < b{. Assume that the interval enclosure of DF(T), denoted by [DF(Vj], is invertible. 
Let Xq € I and define 

N(F, x , 1) = - [DF (I)]- 1 F(x ) + x Q . 
If N(xq,T) C I then there exists a unique point x* € I such that F(x*) = 0. 

5.4.2. Rigorous-interval-based numerical results. For our proof we use a 
normal form of order N — 4 change of coordinates (5.18). At this point we stress 
once again that (j> obtained by (5.18) does not need to perfectly align coordinates. A 
numerically obtained polynomial, provided that it aligns coordinates well enough, is 
sufficient to prove the existence of a center manifold using our method, provided that 
assumptions of Theorem 2.4 can be verified. 

We investigate a set 

= Bf x B r u x Bl 

with 



R = V2 • 155 • 10- 4 « 0.176, 
r = 5- 10~ 4 . 



(5.30) 
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Our choice of R by (5.17a) implies that we consider actions / e [0, 0.0155]. 

We first prove that we have a fixed point in applying Theorem 5.14. We take 
I = nf =1 [-25 • 10~ 5 , 25 • 10~ 5 ] C intD* and compute 

iV(F* 0,1) c 

C [-6.24567e - 05, 6.245664e - 05] x [-6.24434e - 05, 6.24435e - 05] 

x [-6.24908e - 05, 6.24908e - 05] x [-5.33554e - 05, 5.33554e - 05]. 

Clearly 7V(i^,0,I) C I, which establishes that L x is in the interior of . 

Next we verify condition (2.5). We take v = V2 ■ 155 • 10~ 4 - V2 ■ 150 • lO" 4 , 
which is equivalent to the ball B^~ v having actions I € [0,0.015]. We subdivide 
B r u x B r s into 9 pieces and cover dB* by 500 small boxes in E 2 . Taking the 9 • 500 
sets, using (5.29) we obtain a bound on the energy 

H{<j>- l {dB* x B r u x B r B )) > -1.500445781623899. 

Taking same type of subdivisions we then show that 

H{4r l (B*-° x §1 x &D) < -1.500445787982231, 



R r- I |31000j 



which establishes (2.5). 

To compute [L>F*(D )], we cover Bf by 31000 boxes in R 2 , B? C (Jiii 
Taking Ij = l c i x B^ x B r s we compute the bound on [DF^{D^)\ as an interval hull 
of all matrices [DF^(I i )] (Th is means that we take an interval matrix \DF^{D^)\ so 
that C [DF+IDj,)] for i = 1, . . . , 31 000). Each interval matrix [DF^iU)] 

is computed using (5.28). Thus we obtain a bound for [DF^(D^)] (displayed below 
with 3-digit rough accuracy; rounded up to ensure true enclosure) 



( [-0.0336,0.0335] 
[-2.15, -2.03] 
[-0.0783,0.0782] 
\ [-0.0782,0.0783] 



[2.06,2.11] 
[-0.0422,0.0422] 
[-0.0559,0.0566] 
[-0.0559,0.0566] 



[-0.0526,0.0521] 
[-0.0826,0.0827] 
[2.43,2.64] 
-0.0962,0.09741 



[-0.0521,0.0526] \ 
[-0.0825,0.0827] 
[-0.0974,0.0962] 
[-2.64, -2.43] 

(5.31) 



,-forw ,,forw 



We take an = Pv = 2 and a v = Ph = 7 = 1, which clearly satisfy (2.13). In our 
application we deal with a single set N p = Nq = D^, which means that for this set 
p = R. With our choice of parameters condition (2.14) clearly holds. 

Based on (5.31) using (2.9-2.12), (2.15) and (2.16) we compute the constants 
fw , Kj? ack , K^ ack , K^ ack , £„, s s , (5", 8 s needed for the verification of 
assumptions of Theorem 2.4. The computed constants are written out in (8.1) and 
(8.4) in the Appendix. 

Finally, using the boxes I C) j we also compute [ir x ,yF(B^ x {0} x {0})] as the 
interval hull of all [t^F^I^ x {0} x {0})] for i = 1, ... ,31000 (displayed below 
with rough accuracy) 

[K Xty F(B? x {0} x {0})] (5.32) 
= [-0.000954908,0.000819660] x [-0.0009549120,0.000819658], 

from which E u , E s are computed using (2.19) (see (8.1) in the Appendix). For com- 
putation of each \n Xiy F^(l Ci i x {0} x {0})] we in fact need to further subdivide each 
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0.0015 r 




FlG. 5.2. A rough sketch of t^x,y,p x 1 {^t^~ V x x ^s)i w hich gives us an idea of the 
actual size and thickness of the investigated region in which we have proved existence of the center 
manifold. 

box I C[ j into nine parts (this is because E u and E s turn out to be our most sensitive 
estimates). Based on all the computed constants we verify assumptions (2.17-2.21) of 
Theorem 2.4. 

The computer assisted part of the proof has taken 3 hours and 49 minutes of 
computation on a standard laptop (it is possible to conduct much shorter proofs, 
but for less accurate enclosures of the manifold than above). Looking at the con- 
stants (8.1), (8.4) written out in the Appendix it is apparent that assumptions (2.17), 
(2.18) of Theorem 2.4 hold by a large margin. The bottleneck lies in conditions (2.20) 
and (2.21). This follows from the fact that the bounds computed in (5.32) are large 
in comparison to r (see (2.19) which binds the two together). This is because far 
away from the origin the 4-th order normal form no longer gives an accurate enough 
estimate on the position of the manifold, and hence the vector-field in the expan- 
sion/contraction direction becomes noticeably nonzero. A simple remedy would be to 
use a higher order normal form, which would allow for obtaining a tighter enclosure 
and also a larger domain. This would require longer computations and use of more 
capable hardware than a standard laptop. Such computations though can easily be 
performed on clusters. 

Finally let us note that the size of the region in which the manifold is found is not 
negligible. In Figures 5.2 and 5.3 we see our region together with the smaller mass 
(Earth) in the original coordinates of the system. Our set D<p is a four dimensional 
"flattened disc", in Figure 5.2 we can see that the disc is not too thick. On our plot 
the set ttx,y,Px i^ 1 (AjO) nes between the two coloured flat discs (blue disc below, 
and green disc above; in this resolution they practically merge with one another). 

6. Closing remarks, future work. In this paper we have given a method 
for detection and proof of existence of center manifolds in a practical domain of 
the system. We have successfully applied the method to the Restricted Three Body 
Problem. The method is quite general. It can be applied to any system with an 
integral of motion which allows for a computation of a normal form around a fixed 
point. The method also works for arbitrary dimension, which makes it a tool which 
can be applied to a large family of systems. 

The strength of our approach lies in the fact that we can investigate and prove 
existence of manifolds within large domains, and not only locally around a fixed point. 
The weakness so far is that the method only establishes Lipschitz type continuity of the 
manifold. In forthcoming work we plan to remedy this deficiency. In our view, since 
we already have established Lipschitz continuity, similar tools combined with standard 
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Fig. 5.3. A rough sketch of tt x ,y (d'~ 1 {dB^~ v X B r u X B r s ~)\ in blue and 

^x.Y (</> _1 x &u x Sj)) in red. We have proved that the manifold is contained in 

<f>~ 1 (B^~ v X X Bl) and that orbits starting from it never leave </> — ^B,? X X BJ) when going 
forwards or backwards in time. 



cohomology equation arguments can be applied to prove higher order smoothness. 

We would also like to mention that the method allows for rigorous enclosure of the 
associated stable and unstable manifolds through cone conditions used in the proof. 
This means that the presented method can be used as a starting point for computation 
of foliations of stable/unstable manifolds, and next a scattering map associated with 
splitting of separatrices. In our future work we plan to conduct rigorous-computer- 
assisted computations of the scattering map for the RTBP in the spirit of [DMR]. 
Such computations can then be used in the study of structural stability or diffusion. 

7. Acknowledgements. We would like to express our thanks to Daniel Wilczak 
for frequent discussions and his assistance in the implementation of higher order com- 
putations in the CAPD library (http://capd.ii.uj.edu.pl). 

8. Appendix. Here we list the bounds needed for the verification of assumptions 
of Theorem 2.4. Below constants were computed using (5.32), (5.31) combined with 
(2.19), (2.9), (2.10) and (2.12) 

E u = 1.909815022732472, E s = 1.909823931307315, 

5 U = 2.434904529896616, S s = 2.434911565550947, 

e c = 0.09796031906285504, e m = 0.09656707906887786, ' ' 

e u = 0.09737656524499766, e s = 0.09735689577043023. 

Note that [DF't'iD^)] and [w x , y F(B? x {0} x {0})] in (5.31), (5.32) are displayed 
with very rough accuracy. Above numbers follow from their precise version from the 
CAPD software. 

From (2.8) we have obtained the bounds c u ,c s (see (2.11)) using the following 
simple estimates. Our matrix C from (2.8) is of the form (see (5.31)) 



C = 



Si ri 
r2 e 2 



For any matrix C = ( £l Tl | gC and any = (9i, 62) for which \\9\\ — 1, using 



2 2 2 2 
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we have 



e T ce 



(ri+ 7-2)0301 +6101+ 6202 

ki+r_| , . ki+r 2 | 
€ — max h mm £j, max h max E\ 

riEri,r2Gr 2 2 ej£ej,i=l,2 ri6ri,r 2 £r2 2 e;S£i,i=l,2 



(8.2) 



The bound (8.2) is easily computable using interval arithmetic and (5.31) 

c" = 0.08050236551044671, c s = -0.08046115109310353. (8.3) 

Here once again the very rough rounding in (5.31) is evident when compared with 
(8.3). 

Estimates (8.1), (8.3) give us 

K forw = 0.3233133031766185, «J? ack = -0.3232720887592754, 

K forw = 2.289103404031357, ^ ack = 2.191595652437820, (8.4) 

K forw = _2.191592853354867, K*? ack = -2.289115357054329. 
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